Energy dependence of directed flow over a wide range of pseudorapidity in Au+Au 

coflisions at RHIC 
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We report on measurements of directed flow as a function of pseudorapidity in Au+Au collisions 
at energies of ySj^„ = 19.6, 62.4, 130 and 200 GeV as measured by the PHOBOS detector at 
the Relativistic Heavy Ion Collider (RHIC). These results are particularly valuable because of the 
extensive, continuous pseudorapidity coverage of the PHOBOS detector. There is no signiflcant 
indication of structure near midrapidity and the data surprisingly exhibit extended longitudinal 
scaling similar to that seen for elliptic flow and charged particle pseudorapidity density. 



PACS numbers: 25.75.-q 

The study of collective flow in ultrarelativistic nuclear 
collisions provides insight into the equation of state, de- 
gree of thermalization, and the early stages of the hot, 
dense matter created. The elliptic flow parameter, V2, 
has been studied extensively over a wide range of colli- 
sion energies and pseudorapidity The directed flow 
parameter, wi, however, has been studied in less detail at 
RHIC energies [11113. 

The PHOBOS detector is composed of several subsys- 
tems (see Ref. for details). The most important for 
this analysis were the silicon multiplicity array, which 
consists of an octagonal multiplicity detector (OCT), a 
silicon vertex detector (VTX), and three annular ring 
multiplicity detectors (RINGS) located on each side of 
the collision point. PHOBOS has the ability to measure 
nearly all charged particles due to the ~ 47r solid angle 
coverage and to record particles with transverse momenta 
down to about 35 MeV/c (140 MeV/c) for pions (pro- 
tons) at 77 = and 4 MeV/c (10 MeV/c) for -q ~ 4-5. 

This analysis is based on data sets for Au+Au colli- 
sions at y/s^ = 19.6, 62.4^ 130, and 200 GeV, as used 
in the elliptic flow study AH data sets were taken 
with the spectrometer magnetic field off, except for data 
taken at 130 GeV, where field-on data was included to 
maximize statistics. Details on event selection and signal 



processing can be found in 0,1^. The results are shown 
in the most central 40% of the total inelastic cross-section 
for which the trigger system was fully efficient at all four 
energies. More information on triggering and centrality 
determination can be found in Monte Carlo simu- 
lations of the detector performance were based cm the 
Hijing event generator ^ and the GEANT 3.211 [9] sim- 
ulation package, folding in the signal response for scintil- 
lator counters and silicon sensors. 

The analysis is based on the anisotropy of the az- 
imuthal distribution of charged particles detected in the 
silicon pads of the PHOBOS multiplicity array. The anal- 
ysis uses a subevent technique where hits produced in one 
area of the detector are correlated with an event plane 
angle found from hits in another region |lOj |. 

Directed flow is quantified by measuring the first har- 
monic, vi, of the Fourier decomposition of the particle 
azimuthal angle distribution. 
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i-(l + ^2v„cos[n(0-7/^fl,)]), (1) 



Tl=l 



where is the true reaction plane angle defined by the 
impact parameter and beam axis. 

Measuring directed flow using a subevent technique is 
associated with several pitfalls because global momentum 
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FIG. 1: Measured directed flow as a function of rj in Au+Au coUisions at y^sTTiv" ~ 19.6, 62.4, 130, and 200 GeV, using the 
mixed harmonic event plane method (open points) overlayed with the standard symmetric rj subevent method (closed points). 
Note the different vertical axis scales between the upper and lower panels. The centrality ranges shown for both methods are 
those which give good mixed harmonic reaction plane sensitivity. For clarity only the statistical errors are shown. 



conservation can produce non-flow correlations between 
the subevent and the particle under study 11]. This anal- 
ysis circumvents this correlation by using subevent win- 
dows that are symmetric about midrapidity [l^ , thereby 
canceling the back-to-back momentum conservation re- 
coil because each subevent is composed of two equal sec- 
tions in the negative and positive r/ hemispheres. The 
subevent regions used in the event plane calculations 
are located in the OCT (1.5 < \r]\ < 3) and RINGS 
(3 < I77I < 5) subdetectors. The OCT subevent is used 
to flnd vi in the RINGS region {\r]\ > 3). Likewise, the 
RINGS subevent is used to find vi in the OCT (I77I < 3) 
region. For each 77-symmetric subevent window, a reso- 
lution correction is applied that was equal to 



where the N and P labels denote event planes found in 
the negative and positive halves of each subevent win- 
dow. The centrality averaged resolution correction for 
the octagon subevent was 4.1 for 19.6 GeV and 3.5 for 
the other three energies. The centrality averaged ring 
subevent resolution correction was 1.9, 3.3, 4.5, and 3.8 
for the 19.6, 62.4, 130 ,and 200 GeV data sets, respec- 
tively. 

This analysis accepted collisions within ±10 cm of the 
nominal vertex position. In this detector region there 
exist holes in the octagon acceptance to avoid shadowing 



the vertex and spectrometer detectors. The detector was 
symmetrized using the procedure described in 0. 

In addition, weights were applied to the pads in the 
symmetrized detector hit map to correct for phase space 
differences between the detector pads, as well as to ac- 
count for instances where more than one track passes 
through a pad. The weighting procedure is the same as 
applied in the analysis of elliptic flow 0. 

Monte Carlo studies showed a suppression of the re- 
constructed flow signal that is dominated by background 
particles that do not carry flow information, as well as 
the loss of sensitivity due to the hit map synimetrization 
and the occupancy correction algorithm. As in the ellip- 
tic flow analysis [3, this suppression is corrected using 
the ratio of reconstructed to input flow from the simula- 
tion. Typical correction levels were in the 25-30% range 
for the results shown. 

In the subevent method described above, while the flow 
itself is measured using symmetric subevents, the resolu- 
tion correction correlates portions of the detector that lie 
in the forward and backward i] regions. Thus, it is possi- 
ble that a small non-flow correlation due to momentum 
conservation affects the final result through the resolu- 
tion correction. In order to estimate the potential size of 
this and any other non-flow correlations contributing to 
the signal, we also analyzed the data using a mixed har- 
monic event plane analysis [T^ . In the mixed harmonic 
analysis, the reaction plane, tp2, is determined using the 
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FIG. 2: Directed flow of charged particles in Au+Au collisions as a function of -q, averaged over centrality (0-40%), shown 
separately for four beam energies. Note the different vertical axis scales between the upper and lower panels. The boxes 
represent systematic uncertainties at 90% C.L., and (Npart) gives the average number of participants for each data sample. 



elliptic flow information and the directed flow signal per- 
pendicular to ip2 (out-of-plane) is subtracted from that 
which is in the plane of 1^2 (in-plane). Since the true 
directed flow signal is in-plane, the assumption is that 
the directed flow signal out-of-plane is due to non-flow 
correlations. 

Specifically, in our implementation of the mixed har- 
monic analysis, ip2 was found in two subevents from 
—3<r]< —0.1 and 0.1 < 77 < 3 and used along with 
the r/-symmetric eventplanes defined above to find 
Vi{ipi,ip2} as outhned in Two V'l event plane angles 
and two "02 angles were necessary in order to find vi in 
all regions of pseudorapidity such that the particle under 
study did not fall into the regions where either ipi or 'ip2 
event plane angles were determined. 

Fig. n shows the fully corrected signal for the directed 
flow at all energies as a function of pseudorapidity for 
both analysis methods. The la statistical errors are 
shown as solid bars. In both methods, the statistical 
errors exhibit a point-to-point correlation due to shared 
event plane and event plane resolution determination. 
The mixed harmonic method gives results which are con- 



sistent with the symmetric subevent method at 62.4, 130, 
and 200 GeV. At 19.6 GeV the mixed harmonic analysis 
results are in reasonable agreement with the symmetric 
subevent method; however, the analyzing power of the 
mixed harmonic method is diminished at this energy due 
to the weak elliptic fiow signal, as well as a very small 
event sample. 

The agreement between these methods implies that the 
reaction plane determined by elliptic flow is the same as 
that determined by directed flow, within errors. This in 
turn means that the flow and the reaction plane that we 
see in Au-|-Au collisions is dominated by a global flow 
of the particles with minimal effects from "non-flow cor- 
relations" . Furthermore, since the V2 reaction plane is 
dominated by 77 near zero and vi by high jryj, this result 
indicates that the reaction plane orientation is consistent 
over the entire pseudorapidity range. 

Fig. 13 shows the results from the symmetric subevent 
method with the 90% C.L. systematic errors. Several 
aspects of the analysis were studied in order to estab- 
lish the systematic errors. These include hit definition, 
hit merging, subevent definition, knowledge of the beam 
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orbit relative to the detector, dN/drj distribution, hole 
filling procedure, consistency of Vi result when rotated 
by 180 degrees, magnetic field configuration and the sup- 
pression correction determination. The systematic er- 
ror from each source was estimated by varying that spe- 
cific aspect within reasonable limits and quantifying the 
change in the final vi as a function of 77. Also, the dif- 
ference between the results from the symmetric subevent 
method and an odd-order polynomial fit to the mixed 
harmonic method was included in the systematic error. 
The individual contributions were added in quadrature 
to obtain the final systematic errors. 

Historically, vi has been defined to be positive (nega- 
tive) at high positive (negative) rj where spectator matter 
is thought to dominate the signal [T^ . We have preserved 
that convention here, although it is important to note 
that the spectator region falls outside of our acceptance 
at the higher energies. Consequently, the regions of t] 
used to find the direction of V'l have varying spectator 
content as the collision energy increases. Thus, it is nec- 
essary to invert the sign of Vi at 130 and 200 GeV in order 
to preserve the sign convention from the lower energies 
and make a direct comparison of the shapes as a function 
of energy, as shown in Fig. |2 

The results in Fig.[51show the evolution oivi as the col- 
lision energy increases. All four energies exhibit a vi sig- 
nal passing smoothly through zero at 77 = as expected, 
indicating that there are no momentum conservation bi- 
ases in the data. The vi becomes more negative with rj at 
each energy, until a "turnover" point is reached, and the 
vi from both 19.6 and 62.4 GeV becomes positive at very 
high pseudorapidities. This turnover at all energies and 
the large signal seen at high \rj\ for the lower energies are 
features uniquely observed by PHOBOS. These effects 
may be due to protons and nuclear fragments taking over 
from pions as the dominant contributors to the directed 
flow signal at high \ri\. 

The results at 62.4 and 200 GeV are in qualitative 
agreement with results from STAR ^ |^ . Both experi- 
ments show vi for an extended region about midra- 
pidity at 200 GeV, while \vx\ increases as \ri\ increases. 
At 62.4 GeV, PHOBOS observes a turnover of the vi 
signal that occurs at smaller pseudorapidity than what 
is reported in the STAR data. This may indicate that 
vi at high 1 7/ 1 is sensitive to the transverse momentum 
range included in the measurement. Recall that PHO- 
BOS measures protons down to ~ 10 MeV/c while 
STAR has a cutoff at - 150 MeV/c. 

Fig. 131 shows the directed flow where data points from 
the positive and negative 77 regions have been averaged 



together and plotted as a function of 77' = I77I — ybeam- 
Since the directed flow curves are odd functions, the neg- 
ative 77 region was multiplied by -1 before the averaging 
was performed to avoid cancelation. Within the system- 
atic errors (shown in Fig. I^J, it appears that all curves 
scale throughout the entire region of 77' overlap, showing 



0.1 4 r • 19.6 GeV 

Q .|2L ▼ 62.4 GeV 

' : ■ 130 GeV 

^ - A 200 GeV 

0.08- 

1- 0.06 1 
> 



0.04; 
0.02; 

0; 

-0.02^ 

-0.04; 



-5 -4 -3 



-2 -1 



FIG. 3: Directed flow, averaged over centrality (0-40%), as a 
function of 77' = |?7| — yteam for four beam energies. The error 
bars represent the la statistical errors only. 



that, within errors, the directed flow exhibits the longi- 
tudinal scaling behavior already observed in the elliptic 
flow and charged particle multiplicity This con- 
firms and expands on an earlier observation of this scaling 
in the directed flow between RHIC and SPS results 0- 

In summary, the pseudorapidity dependence of di- 
rected flow has been measured for several collision en- 
ergies. At each energy, the vi signal is small at midra- 
pidity and grows with increasing \t]\. At very high I77I, 
a turnover of vi is observed, possibly due to protons 
and nuclear fragments dominating the flow signal in this 
range. When studied as a function of 77', vi appears to 
scale throughout the entire 77' overlap region at all ener- 
gies. 
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